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ABSTRACT 

Nova explosions occur on the white dwarf (WD) component of a Cataclysmic Variable stellar system which 
is accreting matter lost by a companion. A Type la supernova explosion is thought to result when a WD, in 
a similar binary configuration, grows in mass to the Chandrasekhar Limit. Here, we present calculations of 
accretion of Solar matter, at a variety of mass accretion rates, onto hot (2.3 x lO^K), luminous (SOL©), massive 
(L25M0, I.35M0) Carbon-Oxygen WDs. In contrast to our nova simulations where the WD has a low initial 
luminosity and a thermonuclear runaway (TNR) occurs and ejects material, these simulations do not eject 
material (or only a small fraction of the accreted material) and the WD grows in mass. A hydrogen TNR does 
not occur because hydrogen fuses to helium in the surface layers, and we call this process Surface Hydrogen 
Burning (SHB). As the helium layer grows in mass, it gradually fuses either to carbon and oxygen or to more 
massive nuclei depending on the WD mass and mass accretion rate. If such a WD were to explode in a SN 
la event, therefore, it would show neither hydrogen nor helium in its spectrum as is observed. Moreover, the 
luminosities and effective temperatures of our simulations agree with the observations of some of the Super 
Soft X-ray Binary Sources and, therefore, our results strengthen previous speculation that some of them (CAL 
83 and CAL 87 for example) are probably progenitors of SN la explosions. Finally, we have achieved SHB for 
values of the mass accretion rate that almost span the observed values of the Cataclysmic Variables. 
Subject headings: accretion -binaries: close - cataclysmic variables- stars: individual (CAL 83 CAL 87)- 
supernovae 



1. INTRODUCTION 

It is commonly assumed that supernovae of Type la (SN 
la) are the results of thermonuclear runaways (TNR) in the 
cores of carbon-oxygen (CO) white dwarfs (WD) which are 
members of interacting binary systems, and have accreted ma- 
terial from a companion until their masses reach the Chan- 
drasekhar Limit (CL) and a carbon detonation/deflagration 
occurs (Nomoto, Thielemann, & Yokoi 1984; Leibundgut 
2000, 2001; Hillebrandt & Niemeyer 2000). The binary star 
systems that evolve to this explosion, however, have not been 
identified. Given the importance of SN la, both to our under- 
standing of the evolution of the Universe and to the formation 
of iron peak elements in the Galaxy, we must identify the pro- 
genitors of these explosions. 

The binary scenario was originally proposed by Whelan and 
Iben (1973), and as a result virtually every type of close binary 
that contains a WD has been suggested at one time or another 
However, based purely on observations, most of the systems 
that have been proposed cannot be the progenitors since the 
defining characteristic of a SN la outburst is the absence of 
hydrogen or helium in the spectrum (Filippenko 1997). (We 
note, however, that SN 2002ic, a SN la, did show a shell of 
hydrogen in its spectrum [Hamuy et al. 2003].) One of the 
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first suggestions was a classical nova system (CN), but the 
amount of core material ejected during the nova outburst im- 
plies strongly that the WD is decreasing in mass as a result of 
the outburst (MacDonald 1984; Gehrz et al. 1998; Starrfield 
2003). Other suggestions such as Symbiotic Novae (T CrB or 
RS Oph, for example) can probably be ruled out because the 
supernova explosion would take place inside the outer layers 
of a red giant and hydrogen would be present in the spectrum 
(Marietta, Burrows, & Fryxell 2000; Lentz et al. 2002) 

An important development was the discovery of the Super 
Soft X-ray Binary Sources (SSS: Triimperet al. 1991) and the 
suggestions that they were SN la progenitors (van den Heuvel 
et al. 1992 [VDH92]; Branch et al. 1995; Kahabka & van den 
Heuvel 1997 [KVDH97]). They are luminous, L, - lO^^"^** 
erg s"', with effective temperatures ranging from 30 to 50 eV 
or higher (3-7 x lO^K). VDH92 proposed that their proper- 
ties could be explained by nuclear burning occurring in the 
outer layers of a WD at the same rate at which it was being 
accreted (Steady Burning Hypothesis: Paczynski and Zytkow 
[1978)], Sion, Acierno, & Tomczyk [1979]; Iben [1982], and 
Fujimoto [1982a,b]). Therefore, the mass of the WD was in- 
creasing, it could reach the CL, and then explode as a SN la 
(VDH92). Nevertheless, except for the pioneering work of 
Sion et al. (1979), no stellar evolution calculations have been 
done for massive WDs accreting at the high mass accretion 
rates required to test the Steady Burning Hypothesis. Calcula- 
tions were done for lower mass CO WDs by Iben (1982) and 
Sion and Starrfield (1994) but the luminosities and effective 
temperatures of their simulations were too low to agree with 
the observations of SSS such as CAL 83 or CAL 87 (Greiner 
2000), and it was not clear that the WDs that they studied 
would reach the CL. 

In contrast, we find in the calculations reported here that 
the material can burn at mass accretion rates that are far lower 
than the Steady Burning rate (a few x 10~^Moyr"': Kahabka 
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2002). A hydrogen TNR does not occur because the WD is 
sufficiently hot that hydrogen burns in the layers near the sur- 
face as it is accreted. We, therefore, refer to our results as 
Surface Hydrogen Burning (SHB) to prevent any confusion 
with Steady Burning. In the next section, we briefly describe 
the SSS. In the following sections, we describe our evolution- 
ary sequences and end with a Summary. 

2. SUPER SOFT X-RAY BINARIES 

The ROSAT all sky survey (and later pointed observations) 
established the SSS as a class of X-ray emitting systems 
(Triimper et al. 1991). A recent catalog can be found in 
Greiner (2000), and a discussion of their binary properties can 
be found in Cowley et al. (1998). Two systems have orbital 
periods of 1.04 d (CAL 83) and 10.6 h (CAL 87) which are 
longer than those of typical CNe. Unfortunately, there is a 
large range in the luminosities and effective temperatures re- 
ported for a given system (Greiner 2000). Part of the difficulty 
is that, for the systems studied by ROSAT, only a small region 
of the emitted spectrum is visible in the detectors and large 
corrections have to be applied to the observations (KVDH97). 

The tabulated luminosities and effective temperatures of the 
two best studied sources are: (1) CAL 83: L = (0.7-6) x 
lO^^erg s-i and Tes = (5 -6) x lO^K (Parmar et al. 1998). The 
analyses of XMM Reflection Grating Spectrometers spectra 
of CAL 83 imply a temperature ~ 5 x lO^K (Paerels et al. 
2001). (2) CAL 87: L = (3-5) x 10^<^erg s'^ and Teff = 
(6-9) X lO^K (Parmar et al. 1997). The WD masses are, 
as yet, undetermined but the temperatures and luminosities 
imply WD radii. Observed properties for a number of SSS 
(values taken from the above references or Greiner 2000) are 
shown in a theoretical HR diagram (Figure 1) along with the 
results of our calculations (described below). The systems 
shown in this figure are CAL 83, CAL 87, RX J0019.8H-2156 
(RXJ0019), RX J0925.7-4758 (RXJ0925), IE 0035.4-7230 
(1E0035), and RX J0439.8-6809 (RXJ0439). RXJ0019 and 
RXJ0925 are in our Galaxy, CAL 83, CAL 87, and RXJ0439 
are in the LMC, and 1E0035 is in the SMC. 

While CAL 87 appears to be less luminous than CAL 83, it 
is an eclipsing binary with an eclipse depth of 2 mag (Cow- 
ley et al. 1990). An XMM observation suggests that we are 
seeing only the accretion disk in X-rays and the WD must be 
more luminous (Orio et al. 2004). A factor of two increase 
in luminosity would place it on top of our results for 1 .35M0. 
The WD is probably this massive (Cowley et al. 1998). 

3. SURFACE HYDROGEN BURNING SEQUENCES FOR SSS 

We use the one-dimensional, fully implicit, Lagrangian, 
hydrodynamic computer code described in Starrfield et al. 
(1998, 2000) with one major change. The nuclear reaction 
network used in those papers was that of Weiss and Tru- 
ran (1990). In this work, we use the pp-ncnoH-rp network 
of Timmes ^. We also use the equation of state of Timmes 
(Timmes and Arnett 1999; Timmes and Swesty 2000) and the 
OPAL opacities (Rogers, Swenson, & Iglesias 1996). We in- 
clude the accretion energy as described in Shaviv and Star- 
rfield (1988) although it is far smaller than the nuclear energy. 

Our initial model is obtained by evolving a WD through a 
CN outburst. Once all the ejected material is expanding faster 
than the escape speed, has reached radii exceeding lO^^cm, 
and is optically thin; we remove it from the calculations. We 
then take the remnant WD and aUow it to cool to 1O~^L0 
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before restarting accretion. We do this for three nova cycles, 
take the remnant WD from the last cycle, and then begin ac- 
cretion when the WD luminosity has decreased to ~30 L©. 
(This luminosity is that found for VI 974 Cyg three years af- 
ter outburst: Shore et al. 1997). We rezone our initial model 
so there are 95 zones covering the entire WD (our results are 
independent of the number of zones but not the zone mass: 
see below). The initial model has an outer boundary mass of 
1O"22m0, a surface zone mass of 10"^Mq, and a temperature 
exceeding 7 x lO^K. (We report the effects of the initial lu- 
minosity on the evolution in a later paper.) We used this pro- 
cedure to calculate initial models for 1.25Mq, and 1.35Mq 
WDs. The core composition was 50% carbon and 50% oxy- 
gen. Solar composition material was accreted onto these WDs 
at a range of mass accretion rates. The metaUicity of the LMC 
is about one-third that of the Solar neighborhood so we will 
study the effects of metallicity in a future paper. This will be 
important since metallicity clearly affects the evolution of SN 
la (Timmes et al. 2003), and the range of z. over which SN 
la are currently being studied impUes that they have different 
metallicities. 

We begin the evolution with a bare CO core. In Sequence 4, 
for example, we find that with a surface zone temperature of 
7 x lO^K and density of 10"* gm cm"^, the energy generation 
(e„,„ ) rises quickly to 10"* erg gm"'s"^ The outer layers are 
strongly heated by this energy release and, after 15 yr of evo- 
lution, Tsz (the temperature in the surface zone) = 2.1 x lO^K, 
inuc = 3.5 X lO'^erg gm^'s"', and Lshb = 1-3 x 10^*^ erg s"^ 
It is now sufficiently hot in the surface zone so that it takes 
less time than the time step (~ 2 x lO^s) for all the infalUng 
hydrogen to bum to helium in this zone. In addition, some of 
the helium is already burning to carbon in this zone. We des- 
ignate this process Surface Hydrogen Burning (SHB). Similar 
behavior, hydrogen burning out in the surface zone, occurs in 
aU sequences Usted in Table 1. The heUum mass fraction de- 
clines to zero somewhat deeper into the WD and would have 
too small an abundance to be seen in the spectrum if this WD 
were to explode. 

We find SHB at both WD masses but, for brevity, present 
only part of the results for 1.35M0 in Table 1 and Figure 1. 
The results for other masses, and a detailed discussion of all 
our evolutionary sequences, will appear elsewhere (Starrfield 
et al. 2004). The rows in Table 1 are the mass accretion rate 
(Mq yr"'), the length of time in years that we have followed 
the evolution, the amount of mass accreted (Mq), the tem- 
perature (Tsz) and rate of energy generation (csz) in the sur- 
face zone (sz), the luminosity, and the effective temperature 
(both in K and eV) of the evolutionary sequence. They are 
tabulated at the time we end the evolution. We stopped the 
evolution in Sequences 3, 4, and 5 when the central density 
approached values where carbon burning, not included in our 
network, should have begun. Sequence 1 exhibited a TNR in 
the helium layer after ^ 5 x 10"* yr of evolution but insuffi- 
cient energy was released to directly eject any material. This 
sequence achieved a peak temperature in the TNR exceed- 
ing 7 x lO^K and will be redone with a larger network. The 
other evolutionary sequences were hotter and not as degener- 
ate when helium burning began so that no flash occurred. The 
final rows give the composition at the interface between the 
accreted matter and the core matter (CI) near or at the end of 
evolution. 

Our results show that SHB occurs at I.SSMq for values 
of the mass accretion rate typically observed in Cataclysmic 
Variables (CVs). This implies that if a post-CN system, with a 
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massive WD, can arrive at a configuration where mass transfer 
onto the WD begins when the WD is still hot and luminous, 
then the mass of the WD will grow even if it is accreting at 
a rate normally observed only for CVs. The system does not 
have to be accreting at or near to the canonical Steady Burning 
rate for the WD to grow in mass to the CL. 

Our most important result is that in none of our sequences 
is there any hydrogen (and little helium) left in the surface 
mass zones on the WD. Additional sequences done with finer 
mass zoning at the surface (1O~^M0 or 100 times smaller than 
the sequences in this Letter) show that hydrogen and helium 
reach only to a depth of ^ 1O~^M0. This implies that if such a 
structure were to explode in a SN la outburst, then there would 
be insufficient hydrogen or helium to appear in the spectrum 
(Lentz et al. 2002). 

We also find that for sequences 2, 3, and 4, as the mass 
accretion rate increases, the abundance of '^C decreases and 
that of '*'0 increases. In Sequence 1 the helium TNR was 
responsible for fusing the nuclei to the top of the included 
network (mass 20) and Sequence 5 was sufficiently hot to burn 
almost all the nuclei to this mass. Therefore, we predict that 
the C/O ratio in the topmost layers of the ejecta (observed first 
and moving the fastest) will vary between la explosions and, 
when measured, would provide an estimate of the accretion 
rate prior to the la explosion. 

The accretion rate in Sequence 5 is above the assumed 
Steady Burning rate but it still undergoes SHB and does not 
grow rapidly to large radii and shut off accretion. A sequence 
was evolved with an accretion rate of 1.6 x 10"^Mq yr"^ and 
it expanded to large radii after 30 years of evolution. 

Much larger amounts of mass were accreted onto these 
WDs than in our previous studies of accretion onto low lu- 
minosity WDs done for studies of CN outbursts (Starrfield et 
al. 2000, and references therein). In one case (Sequence 5), 
the mass of the WD increased from 1.35M0 to ~1.38M0 and 
the central density grew to 2 x lO'g cm"^. Clearly, these cal- 
culations need to be extended with carbon burning reactions 
included in the nuclear reaction network. 

We compare the results of our calculations to the observa- 
tions in Figure 1 which plots both the Luminosity (Lshb) and 
effective temperature of the sequences together with the ob- 
servations in a theoretical HR-Diagram. The final value for 
each of our sequences is plotted as an asterisk (*) and they 
are connected by a line to guide the eye. The observed ranges 
for the SSS are plotted as boxes. CAL83 and 1E0035.4 fall 
close to the results for L25Mq, while CAL 87 and RXJ0925 
fall close to the results for I.SSMq. We predict that these 
four systems probably contain massive, hot WDs accreting at 
high rates and growing to the CL. While CAL 87 falls be- 
low the Une for 1.35M0, because it is ecUpsing the WD is 
probably more luminous. The other two SSS are more prob- 
lematic. RXJ0019 and RXJ0439 are in a region of the HR 
diagram suggesting that they either contain lower mass WDs 
or are recovering from an extended traverse to cooler effective 
temperatures. Given the uncertainties in the observations and 
the fact that we have only begun this study, the agreement is 
encouraging. 

4. SUMMARY 

We have investigated the effects of accretion at high rates 
onto luminous (3OL0), hot (Teff=2.3 x lO^K), and massive 
(I.25M0 and I.35M0) WDs. We find quiescent burning of 
hydrogen and helium in the surface layers and, thereby, the 
growth of the WD to the Chandrasekhar Limit. Quiescent 
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Fig. 1 . — We compare the luminosities (Lshb) and effective temperatures 
of our sequences for I.SSMq and 1.25Mq to tlie observed X-ray luminosi- 
ties and temperatures for well studied SSS. Because the observers only give 
ranges, we plot their data as boxes with a label indicating the system. The 
final values for each of our sequences are plotted as an asterisk (*) connected 
by a line to guide the eye. The line is labeled by the WD mass. 



burning occurs at a large range of mass accretion rates and 
the luminosities and effective temperatures of our evolution- 
ary sequences agree with those observed for the SSS which 
makes it probable that systems such as CAL 87 and CAL 83 
are the progenitors of SN la. We also find, for all but the 
lowest mass accretion rate, that all the accreted hydrogen and 
heUum is transformed to at least carbon and oxygen in the sur- 
face mass zones. If one of our evolutionary sequences were to 
explode as a SN la, then there would be less than ~ lO'^M© 
of hydrogen and helium in the expanding layers which is un- 
observable (Lentz et al. 2002). We also find that the ratio of 
C/O in the surface zones is a function of mass accretion rate 
and should vary from one SN la explosion to another. 

An exciting implication of this study is that a Cataclysmic 
Variable binary system can produce either a CN or a SSS de- 
pending primarily on the luminosity of the WD when accre- 
tion is initiated after a CN outburst. Therefore, we propose 
first, that a CN and a SSS may be two different phases of the 
evolution of the same binary system and, second, that a CV 
system may evolve from a SSS to a CN or the reverse at dif- 
ferent times of its evolution. 

We also find: 

1) If CAL 83 is a 1.25M0 WD, then it is accreting at a rate 
between 2 x 10"^ Mg yr~i and 8 x 10"^ Mq yr"i. This is less 
than the Steady Burning mass accretion rate. 2) If CAL 87 
is a L35 Mg WD, then it is also accreting at a rate between 
2 x 10"** and 8 x 10"*^ M© yr'^. Since CAL 87 is eclipsing, 
we assume that we are observing only a fraction of the energy 
emitted by the WD. 3) Our initial models were computed by 
accreting onto WDs that had just experienced a CN explosion 
but had not yet cooled to low luminosities. If this evolution 
is realized in nature, then we predict a link between CN and 
SSS which results in SN la explosions. 4) Sequences 3, 4, 
and 5 have accreted a great deal more material than found in 
studies of accretion onto cool WDs. We feel confident that if 
we were to continue the evolution for the necessary time, then 
they would reach the Chandrasekhar Limit. Thus, they satisfy 
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TABLE 1 

1 .35Mo Hot White Dwarf Evolutionary Sequences " 



Sequence 2 3 4 5 



M (M0 yr ) 


1.6 X 10"^ 


1.6 X 10"** 


1.6 X 10"^ 


3.5 X 10"^ 


O A , , 1 7 

O.U X lU 


Tevol (yr) 


8.3 X 10^ 


2.2 X 10^ 


1.3 X 10^ 


3.8 X 10* 


4.8 X 10 


5Macc(M0) 


1.3 X 10"^ 


3.5 X 10"^ 


2.1 X 10"^ 


1.3 X 10~^ 


3.8 X 10 


T,^(10«K)'^ 


114 


177 


319 


347 


516 


esz (10** erg gm^'s"') 


0.15 


1.5 


16.6 


36.8 


91.0 


LsHB (erg s"') 


5.4 X 10" 


5.9 X 10^* 


6.3 X 10" 


1.4 X 10^^ 


3.6 X 10^^ 


Teff(K) 


3.6 X 10^ 


6.6 X 10^ 


1.2 X 10* 


1.4 X 10* 


2.0 X 10* 


Teff(eV) 


32 


57 


107 


125 


175 


'H(CI)'' 


0.0 


0.0 


0.0 


0.0 


0.0 


"HeCCI) 


<0.01 


0.0 


0.0 


0.0 


0.0 


'2C(CI) 


0.11 


0.45 


0.31 


0.22 


0.06 


'^C (CI) 


0.01 


0.05 


0.13 


0.14 


0.12 


"•n (CI) 


0.02 


0.35 


0.13 


0.10 


0.06 


(CI) 


0.01 


0.15 


0.42 


0.49 


0.07 


A>19(CI) 


0.77 


<0.01 


<0.01 


<0.01 


0.60 



"The initial model for all evolutionary sequences had Mwd=1.35Mq, Lwd=30Lq, Tj.ff=2.3 X 10^ K, and Rwd=2391 km 
''Sequence 1 expeiienced a hehum TNR in which the temperature exceeded 7 X lO^K 
'^sz = surface zone 
CI= Composition Interface: All abundances are mass fractions. 



the conditions necessary to be considered as strongly viable 
candidates for the progenitors of SN la explosions 
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